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Abstract 

The unstable manifold of the B = 2 sector of the Skyrme model is constructed 
numerically using the gradient-flow method. Following paths of steepest de- 
scent from the B = 2 hedgehog, we apply a collective coordinate description 
for the motion on the manifold to extract a Hamiltonian, approximately valid 
for the quantum description of the low-energy nucleon-nucleon interaction. 
The resulting potential - obtained in the Born-Oppenheimer approximation 

- is in qualitative agreement with phenomenological potentials. 
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Because of the non-abelian character of QCD, a first-principles derivation of low-energy 
hadronic physics has proven exceedingly difficult. Based on the large- N c limit, where N c is 
the number of colors an interesting non-perturbative approach is offered by effective 
chiral meson models, as witnessed by their great phenomenological success. A popular real- 
ization is the Skyrme model - rather simple in nature yet remarkably realistic JjJ. Its virtue 
lies in the fact that the non-linear classical field equations allow for topological solitonic so- 
lutions, called skyrmions, which carry baryon number, B [|J, thus unifying the mesonic and 
baryonic sector in a common description. Aside from single-baryon properties, the multi- 
baryon problem is of great significance in an attempt to understand complex nuclei within 
this framework. This poses a difficult problem since the starting point is inherently classical 
and 'requantization' is highly non-trivial. Most promising, in this respect, is the B = 2 
sector, especially the description of nucleon-nucleon scattering, which will be the focus of 
the present work. 

The Skyrme model is specified by the Lagrange density 

f2 i 2 f2 

C = ^Tr[d,Ud^} + —TrUrfdpU, tfd v U\? + ~^Tr[U - 1], (1) 

where UeSU(2) denotes the non-linear pion field. The first term is the familiar non-linear 
sigma model, while the second is one of the possible four-derivative terms, derived sys- 
tematically in chiral perturbation theory. A representation of the £7-matrix is of the form 
U = (a + iT7?)/f n , with f being the three Pauli matrices. The parameters of the model 
are the pion decay constant /„- (=93 MeV), the Skyrme parameter e s (=4.76) , related to 
the p7T7r-coupling constant || and the pion mass m n (=138 MeV). The values quoted in 
brackets, which will be used in our calculations, ensure the correct asymptotic behavior of 
the fields, rendering the correct one-pion exchange tail of the nucleon-nucleon potential ||. 

In order to assess the applicability of the Skyrme model beyond the B = 1 sector, where 
one gets reasonable results for the properties of the nucleon [[!]] , a study of the B = 2 sector 
should yield an acceptable description of nucleon-nucleon scattering as well as the deuteron 
bound-state. Starting from classical skyrmion-skyrmion interactions one has to decide on 
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the appropriate field configurations in the B = 2 sector. Since the Skyrme-Lagrangian (1) 
is viewed as a model for large- N c QCD, the semi-classical quantization procedure demands 
that we start from the minimal-energy configurations. As has been argued by Manton 
||, the appropriate procedure to exclude undesirable vibrational modes, is to calculate the 
classical fields on the so-called 'unstable manifold' of B = 2 configurations. In this respect, 
an important role is played by the B = 2 'hedgehog', which is of the form U = e lT7TF ( r \ 
characterized by the 'profile function' F(r) with boundary conditions: F(0) = 2tt and 
F(oo) = 0. This special configuration is a saddle point on the hyperplane of minimal 
energy configurations and can therefore serve as a 'source' for further, less energetic B = 2 
configurations. A complete class of configurations can be obtained numerically via the 
gradient-flow method 0. One evolves the fields along a path of steepest descent, originating 
from the B=2 hedgehog which is initially excited by an unstable mode. The set of all decay 
paths defines a twelve-dimensional unstable manifold, called M 12 , to be identified as the 
manifold of minimal-energy configurations in the B = 2 sector. The unstable modes of the 
B = 2 hedgehog have been reported in |KJ by treating the hedgehog on compactified spatial 



'hypersphere' and have been recalculated in [|T1|] in 'flat' space, using lattice discretization. 
In the latter work the gradient-flow method and its numerical implementation was also 
developed in detail. 



The calculations in [|Tl|] were performed on a three-dimensional grid using a finite dif- 
ferencing scheme for solving the gradient-flow equation. Applying a collective-coordinate 
description to the motion on M 12 , which is always possible locally, the gradient-flow equa- 
tion in collective coordinate space reads as 

The evolution parameter t is a pseudo-time, having a dimension (time) 2 . The matrix Mij(Q) 
is the metric- or mass tensor on M 12 , while V(Q) denotes the interaction potential between 
two skyrmions. Both quantities depend, in general, on three relative collective coordinates 
of the B = 2 system. For well-separated skyrmions, these coordinates can be defined as 



the distance R and the relative isospin orientation C between the skyrmions. Taking the 
separation axis - in the body-fixed frame - along the z-direction, the matrix C can be written 
in terms of two Euler angles (3 and \ as || 

C = cos(^) cos(|) + ir 2 sin(^) + ir 3 cos(^) sin(|). (3) 

In addition to the three relative coordinates, a, B = 2 configuration is characterized by nine 
global collective coordinates, which fix the position of the center of mass and the overall 
orientation in space and isospace. Subsequently, we shall treat these coordinates in the 
adiabatic approximation, since they describe the 'zero-mode motion' of the skyrmions. 

By exciting the B = 2 hedgehog with different linear combinations of the elementary 
unstable modes, one can - in principle - construct the entire unstable manifold numerically. 
Denoting the six unstable modes of the B = 2 hedgehog by 8U^ M ' 1 ' (magnetic mode [fL0| ) 
and 5U( E '^ (electric mode |[10||), where i = x,y, z specifies the symmetry axis of the unstable 
mode (for details see JllJ), a general combination of the unstable modes, in the body-fixed 
frame, can be written as 

5U = (cos ip)5U {M > z) + (simp) (cos6)5U {E ' z) + (sm6)5U ( - E,x A . (4) 



All other combinations of the unstable modes are related to one of the combinations above 
plus a global transformation of the B = 2 hedgehog. Thus, in the body-fixed frame, only the 
two mixing angles ip and 9 are of importance. Together with the pseudo-time parameter t, 
these are three collective coordinates which uniquely specify the field configurations along the 
paths of steepest descent. As discussed above, another complete set of collective coordinates 
is given by the three relative collective coordinates R and C. In numerical calculations the 
latter are easily extractable only for well-separated skyrmions. 

Using the finite-difference method as in [Tl|], we calculate the paths of steepest descent 



from the B = 2 hedgehog numerically. Several paths are shown in Fig. 1 and Fig. 2, where 
the interaction potential along a given path is displayed as a function of the distance R 



between two skyrmions. For the distance coordinate we use a quadrupole definition [JTT 



The computational details can be found in []12~| |. Depending on the initial mixing angles in 



(|j), the excited B = 2 hedgehog, positioned at the maximum of the interaction potential, 
can either split into two distinct B = 1 skyrmions (Fig. 1), or the paths end in the field 
configuration of minimal energy - the well known torus configuration (Fig. 2). The latter is 
classically bound by 120 MeV and is a source for large attraction in the interaction potential 
f|. The path of steepest descent which results from the gradient-flow of two well-separated 
attractive skyrmions {(3 = 7r) is also shown in Fig. 2. The field configurations along this 
path can be described entirely by the relative distance coordinate which, together with the 
nine global collective coordinates, lead to a ten-dimensional submanifold of M 12 , usually 
called M 10 i,|13||[4 . 



Figs. 1 and 2 indicate that the motion on the unstable manifold is generally quite compli- 
cated. The skyrmions move relative to each other in R and also tend to rotate in C whenever 
necessary in order to approach a configuration with less energy. As discussed in [[Uj 



we 



can identify the relative isospin orientation C between skyrmions for three different paths of 
steepest descent on M 12 . Along the two elementary decay paths from the B = 2 hedgehog 
for excitations with the unstable modes superimposed by mixing angles of (ip = 0) and 
(ip = |,0 = 0), the skyrmions have fixed relative isospin orientations C = 1 and C = ir^, 
respectively. Along the gradient-flow path originating from two-well separated attractive 
skyrmions, the B = 2 configurations retain their relative isospin orientation, C = ir 2 . 
Fig. 3 shows these three important paths on the unstable manifold M 12 . Performing a 
finite expansion of the potential energy in terms of Wigner-D-functions, as explained in || , 
the skyrmion-skyrmion potential can be written approximately as 

V(R, (3, x) = V 00 {R) + V 10 (R) cos(/3) + V U (R) cos 2 (^) cos( X ), (5) 

where the radial coefficients Voo, Vio and Vn are deduced from the three potentials shown in 
Fig. 3. We expect that the expansion will hold for medium and large distances. At shorter 
distances, higher-order terms in the expansion, may become important, however. 

Inspection of the element Mrr(R,C) in the reduced mass tensor along the three im- 



portant paths on M , using the method of [11], shows that this quantity depends rather 



strongly on the collective coordinates at shorter distances. In order to keep the reduced 
mass fixed at its asymptotic value fi s = M s /2, where M s = 1463 MeV is the single- skyrmion 
mass in our parameter set, we change the definition of the distance R on M 12 . Different 
choices are possible since, a priori, no distance definition is preferred. The only constraint 
is that all definitions lead to the same result for well separated skyrmions. For the three 
relevant paths the new distance can be easily calculated [12| . The results for the interaction 



potential are also shown in Fig. 3. Due to the redefinition of the distance the potential @ is 
modified at shorter distances. This indicates further that a detailed knowledge of the metric 
structure on M 12 may become of importance for a proper description. 

Since a more thorough calculation of the mass tensor Mij(Q) on M 12 is beyond our 
means at present, we use its large distance form, in order to define a collective coordinate 
representation of the dynamics. The approximate Lagrange function in the body-fixed frame 
then becomes 

L bl = ^R 2 + ^($ 2 + cos 2 (^)x 2 ) - V(R,(3,x). (6) 

The quantity A s = (265MeV)~ 1 is the moment of inertia of a B = 1 skyrmion in our pa- 
rameter set. The quality of the approximation can be tested by solving the corresponding 
gradient-flow equations for different initial values of R, (3 and \- As a result, we find quali- 
tative agreement of the resulting flow curves with paths of steepest descent on the unstable 
manifold. Unfortunately, a precise comparison is not possible at the moment, because of the 
difficulty of a rigorous definition of the relative collective coordinates for all paths on M 12 . 

Treating the nine global collective coordinates in an adiabatic approximation and using 
the large-distance mass tensor on M 12 , we arrive - after canonical quantization of all col- 
lective coordinates - at a simple Hamiltonian for the B = 2 sector. In the cm. system it 
reads 

H = f| + Sr- + 2^ + 3> + ™- + n*. c). (7) 



Here Pr denotes the radial momentum operator, L is the angular momentum operator and 
Si t 2 are the spin operators of the two interacting baryons. As has been done earlier, we 
apply a Born-Oppenheimer approximation to the collective coordinate dynamics, in order to 
deduce a nucleon-nucleon potential from the Hamiltonian above [fig] . Taking into account 
finite- N c effects [|], the reduced Hamiltonian is diagonalized for given R in the space of 
asymptotic nucleon and A-isobar states. This leads to a nucleon-nucleon potential of the 
form 

v}f N \R) = V C °(R) + V^R){t x t 2 ) + V®(R)(oia 2 ) + V^R)^)^) 

+V£(R)S 12 + V±(R)S 12 {t x t 2 ), (8) 

where S X2 denotes the tensor operator. Fig. 4 shows as an example the central part V®(R) 
in comparison to the phenomenologically successful Argonne potential The potential 



deduced from our nuclear Hamiltonian renders almost all features of the realistic potential. 
It has the correct asymptotic strength, reproduces a medium-range attraction and builds up 
a repulsive core at shorter distances. This demonstrates that the Skyrme model can provide 
sizable attraction between nucleons and therefore might serve as a good model for the study 
of bound states. 

In summary, we have presented the first calculation of field configurations on the unstable 
manifold in the B = 2 sector of the Skyrme model. Applying a collective coordinate descrip- 
tion, we were able to define a simple nuclear Hamiltonian, approximately valid at medium 
and large distances. When performing a calculation of the nucleon-nucleon potential within 
the Born-Oppenheimer approximation, good qualitatively agreement with phenomenological 
potentials for the central part has been found. To improve on the results, one has to examine 
the metric structure on the unstable manifold M 12 , in detail. It has been demonstrated that 
such a analysis is numerically possible. 
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FIG. 1. The interaction potential of two skyrmions as a function of the quadrupole distance 
R along paths of steepest descent from the B = 2 hedgehog. Shown are paths resulting from a 
'parallel superposition' of the unstable modes (9 = 0), using different mixing angles if (in degree). 
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FIG. 2. A different set of paths of steepest descent from the B = 2 hedgehog. Shown are paths 
resulting from an 'orthogonal superposition' of the unstable modes (9 = 90 degree), using different 
mixing angles <p (in degree). 



10 




R(fm) 

FIG. 3. The skyrmion-skyrmion potential as a function of the distance R and the relative 
isospin orientation C. The dashed lines give the results using a quadrupole definition for the 
distance while the solid lines denote the results using a redefined distance definition, as explained 
in the text. 
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FIG. 4. The central part of the nucleon-nucleon potential for the Argonne potential (solid line) 
and for the Skyrme model, using the results from the unstable manifold and the Born-Oppenheimer 
approximation (dashed line). 
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